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ABSTRACT 



Aims. We investigate the gas dynamics and the physical properties of photodissociation regions (PDRs) in ICI396A, which is an 

illuminated bright-rimmed globule with internal structures created by young stellar objects. 

Methods. Our mapping observations of the [C ii] emission in 1C1396A with GREAT onboard SOFIA revealed the detailed velocity 

structure of this region. We combined them with observations of the [C i] ^P] - ^Pq and CO(4-3) emissions to study the dynamics of 

the different tracers and physical properties of the PDRs. 

Results. The [Cii] emission generally matches the IRAC 8yum, which traces the polycyclic aromatic hydrocarbon (PAH) emissions. 

The CO(4-3) emission peaks inside the globule, and the [C i] emission is strong in outer regions, following the 8 yum emission to some 

degree, but its peak is different from that of [C ii]. The [C ii] emitting gas shows a clear velocity gradient within the globule, which is 

not significant in the [C i] and CO(4-3) emission. Some clumps that are prominent in [C ii] emission appear to be blown away from the 

rim of the globule. The observed ratios of [C ii]/[C i] and [C ii]/CO(4-3) are compared to the KOSMA-r PDR model, which indicates 

a density of 10"'-10^ cm"^. 

Key words. ISM: lines and bands - photon-dominated region (PDR) - ISM: individual objects: IC1396 



1. Introduction 

IC 1396 is a large Hii region in the Cep OB2 association, ex- 
cited by the 06.5V star HD 206267. Its border includes many 
bright-rimmed molecular clouds (Weikar d et al. 1996 ), among 
which IC 1396 A is the closest globule, 3.7 pc away from 
HD 206267 in the pl ane of the sky assuming a 750 pc dis- 
tance to IC 1396 (Mat thewsiri979h . IC 1396A has a more neg- 
ative velocity ( 8 kms"') compared to HD 206267 and sev- 
eral of the other bright-rimmed clouds, indicatin g that it is lo- 
cated in front of the star and moving toward us (IWeikard et al.l 
Il996l ). IC 1396A thus provides a simple example of photodis- 
sociation regions (PDRs) on a globule, illuminated by a sin- 
gle star from one side. Within the globule a cavity associated 
with the young stellar object (YS O) LkHa 349 is seen in IRAC 
8//m, i^co and ^^COCl-O) maps dNakanoet alj|1989l) . and the 
extinction map obtained from The Two Micron All-Sky Survey, 
2MASS, data (Reach et al. 2009). There are additional cavities 
associated with YSOs, and jReach et al.l |2009) suggested that 
the globule is being reshaped from the inside, with each pro- 
tostar residing in a "compartment" that it has blown out via 
its outflow. Therefore, high spatial resolution, as provided by 
the Germa n REceiver for Astrono my at Terahertz Frequencies 
(GREAlQ iHevminck eFaP 1201 2l) onboard the Stratospheric 
Observatory for Infrared Astronomy (SOFIA; iBecklin & Gehrzl 



' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA / Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung 



120091 : r^ung et al.]l2012h at THz-frequencies is essential for un- 
derstanding the structure and evolution of the IC 1396A globule. 



2. Observation and data reduction 

2.1. [Cii] observations with SOFIA/GREAT 

We mapped the [Cii] emission at 1900.5369 GHz (158//m) with 
GREAT onboard SOFIA during the basic science flights on 22 
and 28 July 2011. On the second flight we observed only the 
[C ii] emission with the L2 channel, whereas on the first flight 
the LI channel was operated in parallel tuned to CO(12-ll) at 
1381.9951050GHz. The observations were made in total-power 
on-the-fly (OTF) mode, with 1 second integration time at each 
dump, and 8"step size. The OFF position is located east and 
outside of the globule at 2lh37'"27^9, 57°31'10" (J 2000). 

Calibration was made by the standard pipeline (iGuan et al.l 
2012). The forward efficiency is 0.95, and the beam efficiency 
is 0.5 4 and 0.51 for LI and L2, respectively (IHevminck et al.l 
I2OI2I) . Because the emission line is sufficiently narrow, only 
a linear baseline was fitted between -30 and -20 kms"', and 
between and 20 kms"'. After excluding data suffering from 
strong standing waves and discarding spectra with excess noise, 
the data were spectrally resampled to a 0.5 kms"' resolution 
and spatially to 25" (0.09 pc at the distance of 750 pc) to obtain 
a better signal-to-noise ratio (S/N). In the following, we use the 
data from the XFFTS-backend, because the othe r backends show 
fully consistent and hence redundant data (see IHevminck et al.l 
I2OI2L for details of the spectrometers). 
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Fig. 1. Integrated (-5 to -15 km s ') intensity map of the [C ii] emission (color), overlayed with the contours of (a) IRAC 8 jjm, (b) 
velocity-integrated [Ci] ^Pi - ^Pq and (c) CO(4-3). The contour spacing is (a) 10 MJy sr"', (b) 2 Kkms"', and (c) 10 Kkms"'. 
The (O,0)-position is 21*'36™50'.7, 57°3riO" (J2000) at LkHa 349. Five asterisks in (b) mark the positions of the spectra shown in 
Fig. [2 HD 206267, the exciting star, is located at (ARA, ADEC)=(17.1', -1.8') and illuminates the globule from the east. 



Fig. \T\ shows the integrated line intensity, summed over the 
velocity range from -5 to -15 kms '. The error in the line- 
integrated intensity, estimated from the rms noise of the base- 
line varies between 2-6 Kkms ' and depends on the number 
of times that the map point was observed. CO(12-l 1) is not de- 
tected, and the upper limit is 4-13 Kkms"' (see the appendix 
for details of the derivation). 



2.2. Complementary observations 

Complementary observations of [Ci] ^Pi - ■'Po-line at 
492.1606510 GHz and CO(4-3) at 461.0407682 GHz were 
performed between Oc tober 2000 and A pril 2001 using the 
CHAMP array receiver (iGustenet al.lll998 l at the CSO. The ar- 
ray with 2x8 pixels was operated in frequency-switching mode 
(throw ±25 MHz for CO and +10 MHz for [Ci], rate 1 Hz). 
CO(4-3) data were recorded on-the-fly, the weaker [C i] in raster 
mode. The auto-correlator back-end array in its high-resolution 
mode provided 2048 channels with 0.15 kms ' spectral resolu- 
tion. Main beams were 15" at 464 GHz and 14.5" at 492 GHz, 
with resp ective main beam efficiencies of 0.52 and 0.51 (see 
iPhilipp et al. 2006 . for more details of the set-up). Both maps 
were convolved to 25" angular resolution. The integrated line 
intensities were obtained in the same manner as for the [C ii] 
emission (Fig.[T]D, c). 



3. Results and discussion 

3. 1. Dynamics of tiie [C n] emitting gas 

Figure [1^ shows an overlay of the [C ii] emission and the IRAC 
8/im emission, which traces polycyclic aromatic hydrocarbons 
(PAHs) in PDRs, showing a good spatial match between the two. 
The [Cii] emission also follows the local cavities originating 
from the YSO outflows mentioned in Sect. [1] One exception is 



the rim at the eastern edge of the globule, where the IRAC 8 fim 
has a strong peak while the [C ii] emission is not prominent. The 
spatial distributions of [C i] and CO(4-3) are different from [C n]. 
CO(4-3) is widely distributed over the globule, having a shallow 
peak in the middle. The [C i] emission is distributed toward the 
southeast and farther out compared to CO(4-3), following the 
ridge-like structure from (ARA, ADEC)~(1', 0') to (0', -2') in 
the IRAC 8yum emission. The [Cii] emission also follows this 
rim, but it is much stronger toward the northern part of the glob- 
ule. 

IWoo tten et al.l (Il983 h showed the intensity distribution and 
the position-velocity diagram of CO(2-l) and '^CO(2-l) along 
the longitude cut that intersects LkHo' 349 (located at offset 
(0,0)). The CO(2-l) distribution is similai- to that of CO(4-3), and 
'^CO(2-l) is similar to [C i] and [C ii] in the sense that the inten- 
sity dr ops around LkHa 349. On the other hand. iNakano et all 
(Il989h showed the cavity near LkHa 349 both in CO(l-O) and 
'^CO(l-O) although it appears as a much steeper hole in the lat- 
ter These differences are compatible with being caused by opac- 
ity: '^CO is more optically thick than '^CO, and '^CO(2-l) and 
(4-3) emissions have a higher optical depth than '^CO(l-O) in 
LTE gas with T > 10 K, although detailed radiative transfer 
modeling is needed to quantify it. 

The channel map of [C n] (Fig. |2|i shows that the gas near 
the rim of the globule is redshifted; going west and toward the 
center of the globule the velocity becomes more negative. With 
the globule in front of the exciting star and moving toward us 
as mentioned in Sect. [1] this velocity gradient implies that the 
western part, which is farther away from the illuminated rim, 
moves faster away from the exciting source. The upper right 
panel of Fig. |2] shows two [Cii] components with significantly 
higher negative velocity (see also the spectra at positions D and 
E in Fig. [3]). In the southern component (position E) the lo- 
cal outflows from the YSOs may also affect the velocity field 
because there are several YSOs in the void-like structure from 



Y. Okada et al.: Dynamics and PDR properties in IC1396A 




2 10-1 
A RA [arcmin] 



2 1 0-1-2 
A RA [arcmin] 



Fig. 2. Channel maps of the [C ii] emission overlayed with the 
contours of IRAC 8 fim emission. 



(ARA, ADEC)~(1.5', 0') to (0.5', -2') dReach et alj|2009l) . and 
the IRAC 8yum emission shows a cometary shape toward the 
southwest. 

On the other hand, the [C i] and CO(4-3) emitting gas does 
not show a significant velocity gradient. Figure |3] shows the 
line profiles at the five positions marked in Fig. [TJ). CO(4- 
3) has a shght blue shift a t position E, which is also indi- 
cated bv lWootten et all (Il983h in the CO(2-l) position-velocity 
diagram along the latitude cut. However, there is a clear 
difference between the line profiles of [C ii] and CO(4-3). 
Considering the spatial distribution discussed above, it seems 
that the gas emitting [Cii] is blown away from the rim. The 
numerical simulations of t he evolution of cometary globules by 
iLefloch & Lazareffl d 19941) suggest that such a cometary regime 
lasts ~ 90% of the cloud's lifetime, and numerous small clumps 
that are ejected and accelerated along the globule appear as 
"fuzz" on the blue side of Lefloch & Lazareff 's model position- 
velocity diagram. The observed structures match this cometary 
regime. 
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Fig. 3. Spectra of [Cn], CO(4-3) and [Ci] at the five positions 
marked in Fig. [1] 



3.2. Comparison with PDR models 

We compared the observed integrated intensity ratios between 
[Cii], [ Ci] ^Pi - ^P p, and CO(4-3) with the KOSMA-t PDR 
model (Rollig et al. 2006) at the five positions shown in Fig. [3] 
To exclude the contribution of the velocity components in [C ii] 
that do not match with the CO and [C i] emission, we fitted the 
[C ii] emission line with the CO(4-3) line profile at each position 
and adopted the integrated intensity of this fit result. At posi- 
tion E, the intensity obtained this way is 57% of the total in- 
tegral of the [Cii]. We used a single-clump model, which has 
as free parameters the mean gas density («), the far-ultraviolet 
(FUV; hv = 6-13.6 eV) flux (x) in units of the Draine field 
(2.7 X 10^ Wm"^), and the mass of the single clump (m) in 
units of the solar mass (Mq). We used two simplified model ap- 
proaches because with three observed lines, there are only two 
independent ratios but three unknowns. In Model 1, we fixed ;i^ 
at the value estimated from the luminosity of the exciting star 
HD 206267, 10^^^ Lq, as an 06.5V star, and the distance to this 
star from each observed position, assuming that 50% of the pho- 
tons have an energy within 6 eV < /zv < 13.6 eV. The calculated 
value of X is listed in Table [T] Then we determined n and m 
by fitting the intensity ratios of [C ii]/[C i] and [C ii]/CO(4-3), as 
listed in Table [T] To be conservative, we estimated the error by 
assuming that each line intensity has a systematic uncertainty 
of 30%. At position B, no combination of n and m can explain 
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Table 1. PDR model results. Values in parentheses show the ranges defined by the uncertainties. 



position 


Model 1 (fixed x) 


Model 2 (fixed log(m) = -1) 


logW 


log(n) 


log(ra) 


log(n) 


logOr) 


A 


2.0 


5.0 (4.5 - 5.4) 


-2.3 (-3.0- +0.8) 


5.2(4.2-6.1) 


2.5(1.4-3.1) 


B 


1.9 


- (4.4-5.1) 


- (0.0- +3.0) 


4.1 (3.8-5.0) 


1.0(0.6-1.8) 


C 


1.9 


4.7(4.3-5.3) 


-0.8 (-3.0- +1.7) 


4.7(4.0-5.7) 


1.9(1.2-2.7) 


D 


1.8 


4.9(4.3-5.1) 


-2.9 (-3.0- +0.2) 


5.0(4.0-5.9) 


2.5(1.3-3.1) 


E 


1.8 


4.9(4.4-5.2) 


-2.0 (-3.0- +0.8) 


5.0(4.0-5.9) 


2.2(1.2-2.8) 



the observed ratio, m is not well constrained because these ra- 
tios are insensitive to the mass; the layers that emit [Cii], [Ci] 
^Pi - ^Pq, and CO(4-3) emissions are not the central core but 
outer layers. With x - 10^ and n = lO^cm"^, Ay at the cen- 
ter is ~ 4 even for a clump with m = lO^^M©. Increasing m 
corresponds to adding layers at deeper Ai,, which does not sig- 
nificantly change the ratio between the [Cii], [Ci] ^Pi - '^Pq, 
and CO(4-3) intensities. Therefore, we fixed log(m) in Model 2 
to a value of -1, i.e. O.IM©, and fitted n and x- A clump with 
m = 0.1 Mo and n = lO^cm"^ has a radius of ~ 0.02 pc, which 
is ~ 5" at the distance of 750 pc, small compared to the beam 
size. The obtained n from both models are consistent within the 
errors, giving lO'^-lO^cm"^ at position B and ~ lO^cm"^ for the 
other positions. The fitted x in Model 2 is consistent with the 
values estimated from the luminosity of the star for Model 1 ex- 
cept for position B. At position B, the very low [C ii]/[C i] ratio 
compared to the other positions is explained by a low UV field 
in Model 2. It implies that the UV radiation is shielded at the rim 
of the globule. 
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4. Summary 

We presented mapping observations of the [Cii] emission at 
1900.5369 GHz (158//m) with GREAT onboard SOFIA in 
IC1396A, which is an illuminated globule with internal struc- 
tures created by embedded YSOs. The [C ii] emission closely 
follows the IRAC 8jum emission, which traces the PAH emis- 
sion. Together with complementary [C i] ^Pi - ^Pq and CO(4- 
3) observations, we investigated the spatial distributions of 
the velocity-resolved emission. [C ii] shows significant velocity 
changes within the globule, while this is not the case for the [C i] 
and CO(4-3) emission. The spatial distribution is also different: 
CO(4-3) has a shallow peak at the center of the globule, the [C ii] 
and [C i] emission is strong at outer region of the globule, but 
their peaks show different locations. The spatial structure and 
velocity distribution is consistent with a scenario in which spa- 
tially unresolved clumps, which emit mainly [Cii], are blown 
away along the globule and are accelerated toward us. A PDR 
model analysis of the observed intensity ratios indicates densi- 
ties of ~ 10^ cm"^ and a UV flux consistent with excitation by 
the external star except for one position, where the UV is indi- 
cated to be low because of the shielding at the rim of the globule. 
IC1396 is selected as one of the key regions to be observed by 
the GREAT consortium because it provides many bright-rimmed 
globules with simple exciting sources, and observations in wider 
regions are expected in future flights. 



Appendix A: Estimate of uncertainties 

We estimate the uncertainty or the upper limit of the integrated 
line intensity from the rms noise (cr) of the baseline as follows. 
cr = 1.1-2.9 K and 0.6-2.1 K in the 25" resolution map for 
[Cii] and CO(12-ll), respectively. At each map position, we 
create an artificial spectrum by generating normally distributed 
(Gaussian) random numbers with a standard deviation of cr. 
Then we integrate this spectrum over the velocity range from 
-5 to -15 kms"', which is the same range to obtain the line 
integrated intensity, and obtain an intensity / (Kkms"'). After 
repeating this trial 20000 times, we obtain a clear Gaussian dis- 
tribution of /, and adopt its standard deviation (0-2) as an un- 
certainty of the integrated line intensity at this map position. If 
3cr2 is higher than the integrated line intensity, we consider it a 
non-detection, and adopt 3cr2 as an upper limit. 
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